perm filename REFLEC.F4[JC,MUS] blob sn#078104 filedate 1973-12-17 generic text, type T, neo UTF8
00100		SUBROUTINE REFLECT(MAX,K)
00200		COMMON FREQ(3,0/50,50),FUNC(50),AMP(50),II(1),IJJ(3000)
00300		DO 100 N=0,MAX*2
00400	100	IF(FREQ(3,N,K).LT.0.0)FREQ(2,N,K)=-FREQ(2,N,K)
00500		DO 201 J=0,(MAX*2)-1
00600		DO 201 I=J+1,MAX*2
00700		IF(FREQ(1,J,K).OR.FREQ(1,I,K).EQ.99999.)GO TO 201
00800		IF(ABS(FREQ(1,J,K)).NE.ABS(FREQ(1,I,K)))GO TO 201
00900		IF(FREQ(1,J,K).GT.FREQ(1,I,K))GO TO 20
01000	C	IF(ABS(FREQ(2,J,K)).LT.0.0001)FREQ(2,J,K)=0.0
01100		FREQ(2,I,K)=FREQ(2,I,K)-FREQ(2,J,K)
01200		FREQ(1,J,K)=99999.
01400		GO TO 201
01500	C20	IF(ABS(FREQ(2,I,K)).LT.0.0001)FREQ(2,I,K)=0.0
01600	20	FREQ(2,J,K)=FREQ(2,J,K)-FREQ(2,I,K)
01700		FREQ(1,I,K)=99999.
01800	201	CONTINUE
01900		RETURN
02000		END